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Abstract 

We critically analyze the claims with regard to the relevance of topological 
phases in the physics of neutrino oscillation made in a recent paper [Phys. 
Rev. D 79, 096013 (2009)] and point out some inappropriate exaggerations 
and misleading statements. We find that the tt phase described in this paper, 
while interesting, is an artefact of two major approximations made in the 
paper. We point out a more robust and more familiar vr phase in the neutrino 
oscillation formulae which can be interpreted as a pure Pancharatnam phase. 
We also make some relevant remarks on the distinction between the geometric 
and the topological phase made in the commented paper. 

PACS numbers: 03.65.Vf, 14.60.Pq 
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In a recent paper pp| it has been claimed that " for the minimal case of 
two flavors and CP conservation, there is a geometric interpretation of the 
neutrino oscillation formulae for the survival and detection probabilities of 
neutrino species". In this paper we first recall the derivation of the main 
result in |T]with a slightly different notation and then show that (i) there is a 
more robust and more familiar tt phase in neutron oscillation formulae than 
the one discussed in pQ which can be seen as a pure Pancharatnam phase, (ii) 
the geometric interpretation and the tt phase discussed in [Tj is an artefact of 
two important approximations made in the paper and thus lacks fundamental 
significance, (iii) the geometric interpretation belongs only to an aspect of 
the neutron oscillation formulae and not to the formulae themselves, (iv) the 
distinction between the topological phase and geometric phase introduced 
in the paper is inappropriate and (v) discuss certain misleading statements 
made in the paper. 

As in [I], let \ua > and > be the two flavor eigenstates represented 
by the two antipodal points lying along the z-axis and let \9,± > he the two 
orthogonal mass eigenstates which lie on the line making an angle 6 with 
respect to the z-axis. Let the initial state be 

\ua>= c+\ei,+ > +c-\ei,- >, (1) 

where the coefficients and i^a- of [Ij have been replaced by c+ and c_, 
the rest of the being the same. This state evolves in time t to the state 

|z/„ >'= e'^+c+|^2,+ > +e'^-c_|^2,- >, (2) 

where and D_ are the dynamical phases acquired during the adiabatic 
evolution of the states \9,± > from |^i, ± > to |02, ± >• 

The survival probability and the transition probability are given 

by, 



Pa = \ < I'aWa >' f = \c+\'^\ < l^a\02, + > |^ + |c_|^| < Ua\92, " > ^ 

+ (c+*c_e*(-^-+^+)< Ua\e2, + >* < Ua\e2, - > +C.C.), (3) 



+ (c+*c_e*(-^-+^+)< iyp\92, + >* < vp\92, - > +c.c.). (4) 
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We first note a rigorous result within the two-flavor model which does not 
depend on any approximations. Since it has been assumed that there is no 
decay, Pa and must add to 1. An inspection of eqns. (3) and (4) makes 
it obvious that this is possible only if the cross terms in the two equations 
add to zero. This implies that the complex numbers Aa — < fa|^2, + >* < 
^a\^2, — > and Ajj — < z/^|^2,+ >* < ^^^1^2, — > should be equal in magni- 
tude and differ in phase by tt, i.e. the phase of the number A^*^!^^ must be 
equal to ±7r. Now 

Ap*Aa = < U0\e2,+ >< Up\92,- >*< Jya\02,+ >* < Jya\02,- > ■ (5) 

Using the fact that < t'/3|^2, — >* =< ^2, —Wis > etc. and rearranging the 
terms we get 

A0*Aa = < I//3|^2, + ><92,+\Ua>< 2^a|^2, - >< ^2, -^/J > ■ (6) 

By Pancharatnam's theorem, the phase of the complex number on the right 
hand side of Eq.(6) is equal to half the solid angle subtended by the closed 
geodesic curve starting at the state \i^i3 >, passing through thestates 16*2, — >, 
l^a >) 1^2,+ > and ending at >, which is a great circle on the Poincare 
sphere. This phase is obviously equal in magnitude to n. This is the well 
known tt phase between the oscillations of intensities of the two different 
flavours which is seen here as a pure Pancharatnam phase of magnitude tt. 
It follows from the condition of unitarity and is an elegant example of inter- 
nal consistency of different principles of physics having apparently different 
origins. Within the two-state model, this phase is independent of any ap- 
proximations. It can be shown easily that this ir phase does not depend 
even on the adiabatic approximation. Let us also note that (a) this phase 
is independent of the phases of the individual states occurring in Eq.(6) and 
(b) this phase can be looked upon as the phase acquired by the state ji^a > 
if it evolved along the closed great circle under the action of a constant uni- 
tary hamiltonian that represents rotation about an axis perpendicular to the 
great circle, i.e. under an SU(2) element that represents a 27r rotation on the 
Poincare sphere about this axis. 

If one stares at Eqs.(3) and (4) for a while it becomes obvious that the 
content of the above result can be exactly simulated by the following polar- 
ization experiment: Let polarization states \02,+ > and \92, — > be incident 
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on the two slits of an interferometer and a polarizer that passes the state 
|z/a > be placed in front of the screen and the position of the fringes noted. 
Let the polarizer now be replaced by one that passes the state {up > and the 
shift in the fringes measured. The above result says that irrespective of the 
incident states the measured phase shift must be equal to n in magnitude. 

The above result has in fact been demonstrated in optical interference 
experiments [2] using the two-state system of light polarization. The results 
shown in Fig. 3 of [2] show that a rotation of a linear polarizer through 90° 
always results in a phase shift of ±7r irrespective of the polarization states of 
the interfering beams. 

We next come to the vr phase discussed by Mehta [Ij. In the cross term 
on the right hand side of Eq. (4), if we substitute for c+ and c_ from Eq.(l) 
we get, after rearranging the terms, the product < z/a|6'i, + > < 6*2, >< 
t//3|6'2, — >< 6^1, — |z/a > multiplying the exponential term. To make this prod- 
uct correspond to evolution of a state along a closed great circle one needs 
two more terms < 02,—\6i,— > and < 6i,+\62,+ > which are missing from 
the product. To compensate for the missing terms, the author first sacrifices 
the arbitrariness of the phases of the individual states in the product and 
then makes the approximation that the hamiltonian is CP non- violating. Let 
us note that the approximation of adiabatic evoution of the states \0i,— > 
and 1^1,+ > to the states |^2, — > and 16*2,+ > has already been made. 
Under these two approximations, the author argues rightly that the phases 
of the missing terms are accounted for exactly if the phases of the pairs of 
states 16*1,— >, 16*2,— > and \0i,+ >, 16*2,+ > in the product are related 
by parallel transport and that the phase of the product is then equivalent 
to that acquired in a unitary evolution along a great circle under a constant 
hamiltonian, i.e. equal in magnitude to tt. If any of the two approximations 
is dropped the result is no more true. In fact the author shows in a separate 
paper [3] that if the adiabatic approximation is retained but the hamiltonian 
is allowed to be CP-violating, the product is no longer equivalent to evolution 
along a closed great circle and the phase of the product is no longer vr but is 
equal to that determined by the solid angle of the distorted curve! The phase 
of magnitude n is thus an artefact of the restriction on the hamiltonian. 

In order to evaluate the significance of the result let us first consider the 
case of evolution in free space or in constant density matter where the states 
16*1 > and |^2 > are the same. In this case the variations of the flavor in- 
tensities are pure sinusoidal oscillations. The sinusoidal oscillation has three 
attributes: (A) amplitue of the oscillation which is determined by the mod- 
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ulus of the cross terms in Eqs. (3) and (4), (B) frequency of the oscillation 
which is determined by the dynamical phase term containing and D_ 
and (C) a constant phase which is equal to in case of the survival proba- 
bility given by Eq. (3) and equal in magnitude to tt in case of the transition 
probability given by Eq. (4). The main results of the paper are concerned 
with only the attribute (C), i.e. the constant phase of the oscillation. Some- 
one who had never heard of the Pancharatnam phase would fix this constant 
phase trivially by requiring that the survival probability and the transi- 
tion probability be equal to 1 and respectively at time t=0. Surely it 
can be seen as a Pancharatnam phase, but the claim on this basis, as in the 
abstract of the paper, that "the neutron oscillation formulae have a geometric 
interpretation" is, in our view, a gross exaggeration. The statement "More 
precisely, the standard result for neutrino oscillations is in fact a realization 
of the Pancharatnam topological phase" on page 9 is also an inappropriate 
exaggeration. It is equivalent to claiming in the first example discussed in 
this comment that unitarity is a realization of the Pancharatnam phase !. In 
our view the frequency of the oscillation determined by the dynamical phase 
and the amplitude of the oscillation determined by the mixing angle are at 
least as important parts of the neutron oscillation formulae as the absolute 
phase of the oscillation. 

When variable matter density is introduced, the dynamical phase is no 
more a linearly varying function of time nor is the amplitude of the cross term 
constant in time. Both these effects are just as important to the variation of 
fiavour intensities along the path as the geometric part of the phase and are 
parts of the neutron oscillation formulae. In fact in the adiabatic limit, since 
the dynamical phase is large compared to the geometric phase by definition, 
the variation of dynamical phase due to the presence of variable matter den- 
sity could easily dominate over the geometric term. Moreover, as discussed 
above, the value vr for the phase found by the author is a consequence of 
the restriction on the hamiltonian and the adiabatic approximation and is , 
therefore, not fundamental. It is just a special value obtained under specified 
constraints. 

A distinction in terminology between the vr phase obtained when CP 
violation is absent and the non-vr phase obtained when it is present has 
been made in the paper, the former being called topological and the latter 
geometric. We find this distinction unnatural and inappropriate since both 
are parts of a single phenomenon and are manifestations of the same basic 
singularity associated with the SU(2) group. A unified description of this 
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singularity has been described in several papers. In the context of adiabatic 
evolution, a three-dimensional generalization of the sign change rule has been 
described in In the context of nonadiabatic evolution an operational 
description of this singularity has been given in |5l El [7] and experimental 
demonstrations with polarization of light have been reported in [61 [21 [8] . We 
briefly recall this work below. 

Consider the evolution of a spin-1/2 state under the action of a hamilto- 
nian which is a function of three parameters x, y and z. First let us consider 
evolution along a closed circuit in a two-dimensional space, i.e. the plane 
?y = 0. Let the hamiltonian be degenerate at an isolated point {xq.Q^zq) in 
this plane. Let the parameters of the hamiltonian be changed so that the 
closed circuit moves from a condition where it does not encircle the point 
(xo,0,zo) to a condition where it does. Then the considerations in [Tj say 
that at the transition point where the boundary of the closed circuit crosses 
the degeneracy, there is a sudden jump of magnitude vr in the phase of the 
state, assuming that the dynamical phase has been subtracted out. We point 
out that at the point of crossing of the degeneracy the adiabatic approxima- 
tion must break down. It was pointed out in |1] that if the same circuit were 
located in the plane y = e, where e is small and the same motion of the circuit 
carried out so that this time there is no actual crossing of the degeneracy, 
there is a measurable phase shift of a magnitude approximately equal to, +% 
whereas if the same operation were carried out in the plane y = —e there is a 
measurable phase shift of a magnitude approximately equal to — vr (only the 
relative sign being important) [9j. It was further shown that if the circuit 
were taken around a closed loop such that it loops the degeneracy, there is 
a measurable phase shift equal to ±2n. A monopole of strength 1/2 located 
at the degeneracy point gives a good unified description of the two effects 
which have been termed differently in [T]. 

To consider the more general nonadiabatic case, an element of SU(2) 
which corresponds to a 27r rotation about any axis on the Poincare sphere is 
represented by the matrix —1. Any state acted upon by this element therefore 
acquires a phase of magnitude vr. This is the well-known phenomenon of "47r 
spinor symmetry". States with different polar angles with respect to the 
rotation axis execute small circles of different diameter but the total phase 
acquired in one full cycle is always of magnitude vr. There is, however, another 
dimension to the problem. A phase shift has a magnitude as well as a sign. If 
the phase of a state evolving under the above hamiltonian were continuously 
monitored with an interferometer with reference to some reference state \R >, 
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the total measured phase shift will be +7r for some states and — tt for others. 
At critical values of the parameters, the phase shift jumps suddenly from 
being +7r to being — tt. This has been verified in polarization experiments 
[6l[2]. This is due to the presence of phase singularities which can be identified 
as follows. Let |-R > stand for the reference state in the interferometer and 
\R > the diametrically opposite state on the sphere which is orthogonal to 
\R >. At points during the evolution where the evolving state is equal to 
\R > the interference pattern vanishes and the phase is undetermined. In the 
vicinity of this point the phase shifts vary sharply. In general, if a state |/ > 
incident on one arm of an interferometer undergoes an SU(2) transformation 
U which is a function of some variable parameters to yield a final state \F > 
which interferes with a state |i? > in the reference arm then a closed cycle of 
the parameters of U around an isolated point at which |F >= \R > yields a 
total phase shift equal to ±2n7r where n is an integer index representing the 
strength of the singularity. Some examples of such phase shifts have been 
demonstrated in interference experiments [8]. This result is conceptually 
simpler than the adiabatic result as it does not depend on subtraction of a 
large dynamical phase from the total phase. Note that in both the above 
discussions the phase shift associated with the singularity is of magnitude 
2nn and not tt. . 

Finally we point out that at the end of page 2 of [1] , the statement "Now 
this open loop (noncyclic) Schrodinger evolution of a quantum state over a 
time T can be closed by a collapse of the time-evolved quantum state at r 
onto the original state at r = by the shortest geodesic curve joining the 
two states in the ray space [10]." is a misrepresentation of historical facts. 
A similar statement made by Samuel on page 960 of [llj is also false. What 
was said in Samuel and Bhandari [TU] was just the opposite. It was stated 
repeatedly on page 2341 of [TU] that the final state can be connected to the 
initial state by any geodesic arc. After a careful discussion with the first 
author of [10] the author of [1] has stated in [12j that the word "any" in "any 
geodesic arc" refers to any of the gauge copies of the geodesic in the N - space. 
Our response to this is that the expression "any geodesic arc" includes all 
geodesic arcs in N, i.e. those that project down to the shorter geodesic as 
well as those that project down to the longer geodesic in the ray space. The 
footnote from [10] cited in [12] is merely a description of a property of the 
shorter geodesic and does not constitute a restriction on the definition of the 
noncyclic geometric phase. In fact the footnote does not form part of the 
discussion of the noncyclic geometric phase on p. 2341 where the expression 
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"any geodesic curve" occurs again. In response to the following comment 
made in [12] : "the author has been caught in the unfortunate position of 
having been scooped by himself, a statement we do not understand, we 
reiterate that the first clear statement that the noncyclic geometric phase 
should be defined as half the solid angle of the area obtained by closing the 
open curve in the ray space with the shortest geodesic arc connecting the 
final state to the initial state was made in [5], i.e. [R. Bhandari, Phys. Lett. 
A 157, 221 (1991)] and not in [10], i.e. [J. Samuel and R. Bhandari, Phys. 
Rev. Lett. 60, 2339 (1988)]. This led to the prediction of observable ±7r 
phase jumps in two-state systems which were later verified in interference 
experiments [HI [2] and to the nonmodular view of the topological phase, thus 
constituting a conceptual advance in the subject. The restriction to the 
shorter geodesic is thus much more than a footnote. It may also be pointed 
out that the definition of the noncyclic geometric phase proposed in ^ as 
the difference of the total Pancharatnam phase of the evolving state and the 
dynamical phase as defined by Aharonov and Anandan does not depend on 
a geodesic rule and is thus particularly useful for systems with more than 
two states where the geometry of the ray space can not be easily visualized. 
An extension of the definition to the case of an arbitrary reference state has 
been proposed in [7]. To end this discussion we note that the fact that the 
shortest geodesic rule was not used in [lOj was also pointed out in [1]. 

To sum up, the net content of [T] would be precisely summarized if the ab- 
stract of the paper read: "We show that, under the adiabatic approximation, 
the phase appearing in the neutrino oscillation formulae has a geometric 
contribution which, under the constraint of CP non-violation, is equal in 
magnitude to tt ". Considering that the phase in a quantum evolution in 
general has a geometric part, this is not very significant. 
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